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Abstract 

We describe the relationship between time-ordered and retarded response functions 
in a plasma. We obtain an expression, including the proper ze-prescription, for the 
induced current due to hard thermal loops in a non-Abelian theory, thus giving the 
non-Abelian generalization of the Kubo formula. The result is closely related to the 
eikonal for a Chern-Simons theory and is relevant for a gauge-invariant description of 
Landau damping in the quark-gluon plasma at high temperature. 
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High Temperature Response Functions and the Non-AbeUan Kubo Formula 

R. Jackiw and V. P. Nair 



1 Introduction 

The behavior of electromagnetic fields in a plasma of charged particles is described by the 
polarization tensor Il^'^{x,y), which is the two-point current correlation function; pertur- 
batively this is a one-charged-particle-loop diagram with two external photon lines. The 
real part of this tensor describes phenomena such as Debye screening and propagation of 
plasma waves; the imaginary part describes the damping of fields in the plasma (Landau 
damping) 0]. If we integrate out the charged fields in a functional integral for the theory, 
the polarization tensor naturally emerges as the thermal average of the time-ordered product 
of two currents. However, there are situations where the response of the plasma to the elec- 
tromagnetic field is described as the average of the retarded commutator of currents. For the 
real part of the response function, the retarded commutator and the time-ordered product 
coincide. We calculate the imaginary part for both and discuss their relationship in terms 
of spectral representations [Q]. 

In a non-Abelian plasma, such as the one arising perhaps from quarks and gluons, the 
response of the plasma to fields is more complicated. At high temperature, a perturbative 
approach can be used. Nevertheless, it is necessary to go beyond the two-point function; one 
must include at least part of the contribution from all higher point functions for reasons of 
gauge invariance. There is also the need for a resummation of perturbation theory [^. One 
must first sum up 'hard thermal loops' to define new effective propagators and vertices. Hard 
thermal loops are thermal one-loop diagrams in which the external momenta are relatively 
soft (~ gT, where g is the coupling constant and T is the temperature, taken to be high) 
while the loop momentum is relatively hard (~ T). One must reorganize perturbation theory 
in terms of effective vertices and propagators defined by these hard thermal loops, before 
integration on small values gT) of loop momenta can be carried out [Q. This is neces- 
sary so that all contributions to a given order in the coupling constant can be consistently 
included. 

Recently it was shown [§] that the generating functional for hard thermal loops, or 
equivalently the effective action which describes Debye screening, plasma waves, etc., is 
given in terms of the eikonal for a Chern-Simons theory, which had been constructed earlier 
in a non-thermal context [Q. We now describe how this result is extended to include some 
of the decay and damping effects. Our discussion is still at the level of hard thermal loops; 
we do not address effects that arise by incorporating effective propagators and vertices in 
subsequent soft momentum integrations. Our result can be considered as a non-Abelian 
generalization of the Kubo formula [0]. 

In Section II we present the time-ordered and retarded polarizations for the Abelian 
plasma at high temperature. The physical significance of the imaginary part is explained 
in Section HI. The induced current due to hard thermal loops in a non-Abelian theory is 
constructed in Section IV. 
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2 Response Functions for the Abelian Plasma 

We consider QED with the fermion current = and interaction Lagrangian Ci^t = 

—A^J^. The equations of motion for the gauge field have J'^ as a source term. The 
calculation of the current is thus a suitable way of studying response functions. The 
scattering operator S is given by Texp[— i / ct^xA^j,]^] where = eipij^ilji is the current in 
the interaction picture, denoted by the subscript /. We can write in terms of S as 

The equation of motion for the field A^ is 

^^F'^^^{x)=^S-'-^ (2.2) 
oA^{x) 

[While no gauge choice is specified in (|2.2| ), it can be fixed for example by adding —^{d-AY 
to the Lagrange density, whereupon the left side would become OA^ (Feynman gauge).] We 
can use (^3) ^'^ obtain J'^ as a power series in A^^. In particular, a Taylor expansion of (|2.1| ) 
to linear order in A gives 



r{x) = fix) -ijd% e{x' - y')[f{x),3''{y)]A,{y) (2.3) 
Using this in ( |2.2| ) and taking a thermal average we can write 

d,F''^{x) = J d'yU^j,''{x,y)A,{y) (2.4) 

where 

n^'^(x,2/) = -^5(a;°-2/°) mx),f{y)]) (2.5) 

The angular brackets denote thermal averaging, with the unperturbed density matrix e~^°^'^, 
so that {jf^) vanishes. (The Boltzmann constant is set to unity.) Thus the response function 
(p.5|), viz. the average of the retarded commutator [or equation (p.4|)], is appropriate to the 
situation where we perturb the plasma by the field and ask how the field evolves. Equation 
( p.5| ) is the Kubo formula. 

The effective action T{A), which is the sum of fermion loops, is given by 



ir{A) ^ ^J.^-^JA,j^^ (2.6) 

sr{A) 



The current that is here naturally defined as J^{x) = — is given to linear order in Afj_ 



by 

J'^{x)= j d'yIi^^{x,y)A,{y) (2.7) 

where 

Ii^\x,y) = -i{Tf{x)3^{y)) (2.8) 

Clearly, the two response functions do not coincide. One difference is immediately evident: 
The time ordered product is symmetric in its labels (/i, z^)(x, y) [as it must be since Iljf arises 
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= / cosp{x - y) 



from the generating functional T{A) ~ 1 — d^xd'^yAf^{x)Ilqy{x,y)A^{y) + ■ ■ ■ ] while the 
retarded commutator does not enjoy any such symmetry. 

To evaluate ( |2.5| ) and (|2.8D, we make use of the free- field commutator. 

(2.9) 

{Po = v^P^i since we are interested in the high temperature contribution, fermion masses are 
negligible by comparison.) The averages needed in the computation can be evaluated using 
the following expression for the propagator. 

S{x,y) = {Tijj{x)i){y)) 

^ ^-^ [e{x^ - y^) (a,e--(^'-) 

- - x°) (/3pe-*P(""^) + ape-'P^''-y^ 
where 

CXp 1 Pp 

(The bar refers to antifermion distributions.) In evaluating the two-point functions, we shall 
use (|2.9|), ( |2.11| ) and carry out the time-integrations first, introducing as needed convergence 
factors e^'^y , e small and positive. We then get energy-denominators with ze-terms. Further 
simplification and extraction of the imaginary part can be easily done at this stage. We find 
for the polarization tensor, 

W{x,y)= I -^^^e-'^^''-y^W{k) (2.12a) 
J (27r) 



where 



d^q 1 1 
(2^ 2^ 2^ 

+ T^'^(p,g')( 
+ T^''{p',q) ' 
+ T'^''ip',q') 



T^'ip.q) 



,Po -qo-ko-ie - ?o - + iv , 

PpPq \ 



^Po + qo - ko - Po + qo - ko + irj^ 

OipOtq I3p(3q 

,Po + qo + ko-ir] Po + qo + ko + ie ^ 

( QpPq j3pOiq 



,Po-qo + ko-ie Po- qo + ko + irj^ 



(2.12b) 



T>^''{p,q)= tr (Yl-pYl-q) 



(2.13) 



q'' = {q',-ci), 



p 



iPl 



q 



|q| 



and p = q + k. 



Further, rj = e for T products, i.e. for 11 j.'^, and rj = —e for the retarded function 11^'^. 
Notice that the retarded function can be obtained by continuing the real part by the rule 



4 



High Temperature Response Functions 



ko —>■ ko + ie. Formula ( p.l2b| ) agrees with Ref. where time-ordered products were used. 
(We have changed the sign of k in our definition of the Fourier transform relative to Ref. 0.) 

The real part of 11^'^ (A;) is the same for both the time-ordered and the retarded functions 
and has long been familiar [jl| . Here we concentrate on the imaginary part which is obtained 
from ( |2.12bD with the help of 



z — 



P- + m5{z) 



(2.14) 



Consider first the contribution due to the T^^{p,q') and T^^[p',q) terms in ( ^.12b ). The 
imaginary parts carry the (5-functions, 5{pQ + go i ^o)- It is easily seen that these are 
subdominant at high temperature, for k small compared to T. [Their contributions are 
0{T) or less.] The dominant contributions to the imaginary part of 11^'^ come from the 
T^''{p,q), T^'''{p\q') terms. For H^'', we find from (gA2b| ) 



Im = TT 



d^q 1 1 
(2^ 2^ 2^ 



T''''(p, q) {rig - Up) 6{po - qo - ko) 
-T'^''{p',q') {ng-np)6{po-qo + ko) 



(2.15) 
(2.16) 



For Hj^ , we find 



Im Ui^" 



Im n^'' + 2n 



d^q 1 1 
(2^ 2^ 2^ 



T^'^ip, q)npil - ng)5{po - qo - ko) (2.17) 
-I- T'^'ip', q')ngil - np)5ipo - qo + ko) 



In the high temperature limit, |k| is small compared to |p|, |q|, and there are simplifications 
i, 

T^''{p,q) c^.SqlQ^Q'' 
T^'ip^q') - ^qlQ'^Q'" 
Po - qo - ko ^ -k ■ Q 
Po - qo + ko ^ k ■ Q' 

with Q'^, Q'^ being the light-like four- vectors 

g^ = (l,g) , Q'^ = (l,-g) 

Equation ( |2.1(j| ) simplifies as 

k^pi^^ roc 



(2.18) 



(2.19) 



Im IId'^ ~ " „ / dq q{nq + fig 

<j 



koT' 
12 



where 



P"" = JdQ 6{k ■ Q)Q''Q' 



(2.20) 

(2.21) 
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and we take n = n. The integration in ( p.21| ) is over the orientations of the unit vector 
Q = g. Note that from its definition, P^'^ is transverse and traceless, while the (^-function 
enforces k to be space-hke. Exphcit evaluation gives 



6n 



where 



and 



^2 



-'2 



3 ^ 2 V 



A;2 



= 



(2.22) 



(2.23a) 



(2.23b) 



Simplifying (2.16) similarly, we find 



im ~ im n^'' + 



p^J■I^ POO 



roc 

/ dqq^nq{l-nq) 
Jo 



p3 

Im n^^ + — P^" 

^ 6 



(2.24) 



This relationship between IIj, and 11^ can be understood in the following way. 11^ 



being retarded, obeys a spectral representation of the form |0 



■'■■'■sub ^ 



dk, 



IJ 



k'o — ko — ie 



(2.25) 



for some spectral function p^^{k). H'^^^ is a 'subtraction term' that can arise in the real part 



of n^'^. For Ili^ , we then have 



■'■■'■sub ^ 



dk, 



iJ 



where 



k'o — ko — ie 
1 



+ 2mfiko)p^''{ko,k) 



(2.26) 



(2.27) 



The bosonic distribution function fiko) appears because Iljf is ultimately part of the bosonic 
{i.e. photon) propagator, and also because it is given by the thermal average of the T-product 
of two bosonic operators: the two currents and j'^. From (|2.25|) and (|2.26|) we find 



Im = np'^'ik) 

Im n^" 



Im n^^ + 27r/(A:o)p^'^(A;) 
71 coth ^p'"'ik) 



(2.28) 



6 



High Temperature Response Functions 



The essence of our results ( |2.20| ) and ( |2.24| ) is that the high-temperature spectral function is 

k 

pM^(A;) ~ -^P"'' (2.29) 



12 



and the difference in the high-temperature behavior between Im U'^ and Im Ut^" [O(T^) 
vs. 0{T^)] is attributed to the presence in the latter of 27i f (ko) p^'^ , which according to ( p.27| ) 
and (ICTI) tends to ^P^"'. 

The spectral function ( p.29|) also determines the high-temperature behavior for the (com- 
mon) real part of 11^'^ and Hj^, apart from possible subtraction terms. Inserting ( ^^.29] ) into 
the dispersion formula ( |2.25|) or (|2.26|) shows that the well-known high-T asymptote for n- 

of ( p:m| ) 



Re n 



1 + 



2|k| 



In 



ko - 


|k| 


ko + 


|k| 



1 ptj.u 

3 ^ 



-pr 



(2.30) 



is reproduced with an appropriate subtraction term. 

Our result (|2.20|) , (|2.3CI| ) for the retarded function 11^'^ agrees with various previous cal- 
culations 1^. It is noteworthy that these early calculations in the Soviet literature, based 
on the Boltzmann and Vlasov equations of kinetic theory, are here regained in quantum field 
theory at one-loop order. 

Another correlation function that is frequently considered is the imaginary time one. It 
too is given by a dispersive integral @]. 



-^imaginary 
time 



■'■■'-sub ^ 



27rmT 



dk 



kn Un 



(2.31) 



Because the "external energy" a;„ is temperature dependent in imaginary time, it makes 
sense to speak of high temperature behavior only for the n = mode, (effectively reducing 
dimensionality to three) where the spectral function enforces an O(r^) large-T behavior. 

The foregoing discussion emphasizes that finite-temperature field theory can be described 
by different correlation functions, with differing large-T behavior. In particular the retarded 
commutator is what is relevant for the operator equations of motion. Landau damping, 
etc. Because of the asymmetry in its time-arguments, the retarded commutator cannot 
be obtained by varying an effective action, which necessarily involves the symmetric, time- 
ordered products. Correspondingly, only the retarded commutators behave as for large T, 
while the time ordered products and the real-time effective action have a more complicated 
large-T behavior. Finally, the imaginary time correlations again behave as T^ for large T, 
provided the "external energies" are restricted to the tt, = mode, or are continued away from 
the temperature-dependent values 2'KinT. Such analytic continuation from the imaginary 
time expression will produce the retarded functions in real time, and not the time-ordered 
ones, which possess non-analytic [~ f{ko)p'^'^] contributions. 

In spite of this variety, a universally true statement can be made about the high- 
temperature asymptote of the spectral function p'^'^, which determines quantities of physical 
interest through the various dispersive representations. In the present context, such a result 
is given in Eq. ( |2.29D . 
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3 Physical Significance of Im 



The imaginary part of n^*^ describes Landau damping, which can occur for fields with space- 
hke momenta [0. Exphcitly we consider 



(3.1a) 



with 



A^{k) = 6{ko - u;{k))A,{k) + 6{ko + a;(k))At (fc) (3.1b) 

and u"^ < k^. The amphtudes v4^(/c), Ajj^{k), respectively of the positive and negative 
frequency terms, correspond to absorption and emission processes. The decay of the field 
in the plasma arises from absorption by fermions and antifermions. The amplitude for 
absorption by fermions is given, to the lowest order in coupling constant, by 



A = lUpij ■ A)ug{2iTYS^''\p -q~ k)^nq{l-np) 



(3.2) 



where Uq are wave functions for outgoing and incoming fermions, respectively. The factor 



Uqil — Up) arises since the initial fermion is chosen from a state of occupation number Ug 
and the final fermion is scattered into a state of occupation number rip. Single quantum 
absorption as in ( |3.2| ) is kinematically allowed for spacelike momenta; it is in fact the inverse 
of Cerenkov radiation. One can also have creation of the mode [u, k) by (Cerenkov) radiation 
from fermions, given by a formula like ( |3.2| ) with p ^ q and A A"^ . There are similar 
contributions from antifermions. For the net absorption probability per unit spacetime 
volume, denoted by 7, we then find, with summation over all fermion states. 



7 



^ — — AUk) 
(27r)3 2po 2go ' 



T^^iP, q) K(1 - np) - np{l - rig)] 2it6{po - Qq - ko 



+ T'^''{p',q') [np{l-fig)-fig{l 



rip)] 27t5{po -Qo + ko 



where po = |p|, qo = |q|. Comparing this with ( p.l2b|) or (|2.16| ) we see that 

7 = 24(fc)[Im U'^^iko)]Aik) 
We can parametrize the field A^ as 

kko 



A 



k2 



■Ao + At 



(3.3) 



A^k) 



(3.4) 



(3.5) 



where k ■ At = 0. gives rise to the transverse electric field and the magnetic field; 
(f) = (1 — p-j gives the longitudinal component of the electric field via = — V(^. In 
terms of the parametrization (pj.5D we can write 



7 = r- 



,71UJ 



-0V 



6 



(3.6) 



7 is positive, as expected for net absorption or damping. 
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4 Induced Current in a Non-Abelian Plasma 

We now consider the non-Abelian plasma and the induced current due to the hard thermal 
loops. 

As far as the two-point function is concerned, there is no significant difference between 
the Abelian case of electrodynamics and QCD. For the latter, the contribution of Np flavors 
of quarks at high temperature is 

Im n'^"''^' ~ (5""^ (^^P^^^ . (4.1) 

The real part of the high-T two-point function is given by |^ 



Re W"'"'' ~ S"''-^ 
24n 



k-Q 



(4.2) 



The singularity when k ■ Q vanishes is defined as a principal value, and then the integral 
reproduces ( |2.30| ), apart from the group factors. 

Notice that ( [l.l] ) and, with appropriate modifications, ( |2.20D are obtained by continuing 



( [4.2| ) by ko — >■ k^ + ie. We expect this to be true in general, including the gluon contributions. 
Thus we expect 

Im n'^"''^' - 6"' (^N + ^) ^^"'^ (4-3) 

for SU {N) gauge theory with Np flavors of fermions in the fundamental representation. 

We now consider the higher point functions in QCD. The evolution of fields in the plasma 
is still described by an equation similar to ( |2.2| ), viz. in the Feynman gauge 

nA^,{x) = tS-'-^^ (4.4) 

This is still the operator equation of motion, reexpressed using the scattering operator S. 
The operator S is now more complicated; correspondingly the current defined by ( |4.4| ) is a 
more involved expression, including ghost and gluon terms. Nevertheless, when we expand 



( [4. 41 ) in powers of A^, we get only retarded terms, since S ^jjSm does not involve fields in 



the future of ^^(x). In fact, writing U{x^,y^) = Texp{i JyO (i^xAnt), -S" = ?7(oo,— oo), we 
find 

Jai^)='S-'-^ = U{-oo,x')j!i{x)U{x\-oo) (4.5) 



SA-ix) 



ja{x). When we differentiate ([4.5|) with respect to the potentials, we 



where - /f^T\ 

get products of currents, and 'contact' terms like ^^r^- Such terms eventually give rise to 

expressions like f'^^AlF^'' needed to write (U) with a gauge covariant derivative D,yF'^^, 
along with ghost and gauge-fixing terms. [Strictly speaking, we need suitable T*-products 
in the definition of 5* to make the contact terms together with the commutator terms come 
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out covariant 0. Our conclusion will not be affected by this needed qualification, since we 
use only the retardation property.] Ignoring contact terms for the moment, we find 

^"^'^""^ - i-^r E 0{x' - y',)e{yl - yl) ■ ■ ■ - y^) 



6Al\{y^)---5All{yn) 

X [■■■mx),3ll{y{)l3ll{y2)]..] (4.6) 

where the summation is over all permutations of {vi hi yi, 1/2 62 y2, • • •)• The right hand side 
of ( [4.6| ) is the definition of the multiple retarded commutator. Using this, we arrive at 



00 „ 

^a^(^) = faix) + T^i-^r / . . . d'yneix' - 2/?)0(i/° - y°) " " " 0{yU - yl) 

n=l •' 

X [. . . Mix), 3ll{yi)l 3?M)] ■ ■ • ■■KliVn) (4.7) 



We must take thermal averages on the right hand side of (4^). If we write 

= / etc., (4.^ 



the retardation property of (4/7) indicates that the integration over the time-components of 
?/j's requires — > + ie for convergence. This will generally be true, since it follows from 
the retardation property and does not depend specifically on the expansion ([4.7|) . We can 
thus make a rule for the n-point functions: Calculate the thermal average of the current, first 
ignoring the ze's in the propagators or energy-denominators. Then for each momentum ki at 
each of the yiS (but not at the unintegrated point x), make the replacement k^ ^ k^ + ie. 
The averaged equation of motion 

^Alix) = (j;(x)) (4.9) 

gives the evolution of fields in the plasma, {J^{x)) being calculated as above. 

We can now apply this argument to the induced current due to hard thermal loops in 
QCD. The hard thermal loop contribution, found in Ref. 0, is determined in the following 
fashion. We begin with the functional I{A) defined on Euclidean 2-space [complex coordi- 
nates z and z] and depending on the Lie-algebra valued vector potential A, with a single 
(spatial) component UTOj . 

_ j ^ (~^)" f ^^^1 ^^^n tr A{Zi,Zi) . . . A{Zn, Zn) 
Z-ij — Z^ Zj 

This is the Chern-Simons eikonal, i.e. it is the exponent in a WKB wave functional. In a 
Chern-Simons gauge theory on (2-|-l) dimensional space-time, the Schrodinger-picture state, 
which is defined on 2-space (at fixed time) and satisfies the Chern-Simons Gauss law, has 
the form e^^^\ with A being one of the two spatial gauge potentials of the theory and 
the other being conjugate to A. Owing to the simple symplectic structure of this theory, 
the eikonal/WKB approximation is exact, and I{A) coincides with the eikonal. /(v4) is 
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also recognized as the Polyakov-Wiegman determinant for two-dimensional, single chirality 
fermions 0]. The structure ( 4.10 ) is relevant to high-T field theory owing to the fact that 
its gauge transformation properties are closely related to those of the generating functional 
for hard thermal loops 0, and the latter is completely determined by its response to gauge 
transformations. 

To present the thermal generating functional ( |4.10|) must be elaborated upon and con- 
tinued from the Euclidean space {z, z) to Minkowski space-time. This is accomplished by 
replacing A in ( ^.10| ) with = ^A- Q, while z is replaced by x ■ Q', z hj x ■ Q. Also all the 



A's ^ ^A ■ Q depend on four variables: z 



X 



Q', z 



X 



Q and x^, which is orthogonal 



to Q and Q' . All the A's carry the same value of x^, which therefore is a fixed parameter in 
the functional ( |4.10| ), and is subjected to a single integration: I{A) J (Pyi±I{A) = X{A) 
The hard thermal loop generating functional is now given by 



2 J Qn 



J dil J (fx tr A+A^ + mT{A+) + mT'{A_ 



(4.11) 



where 1'{A_) is obtained from X(v4_|_) by interchanging Q and Q', and A_ = ^A ■ Q'. 

The integrals in Minkowski space-time are now singular: when transformed to momentum 
space they contain denominators involving k ■ Q and k ■ Q', which can vanish. If these zeroes 
are ignored and the integrals are evaluated formally, one regains the hard thermal loop 
contributions, but they do not have the appropriate analyticity properties, and cannot be 
interpreted as time-ordered products, retarded products, etc. 

Here we give the prescription which results in an evaluation of the retarded current, 
relevant to the non-Abelian gauge theory. 

The contribution of X(v4_|_) in ( 4.11 ) to the current is given in Euclidean space by fU^ 



61(1 A ■ Q) 
SA-ix) 



N + 



n+l 



n=l 



(^) A+(xi)---A+(x„) 



TT 



{z - Zi)zi2 ■ 



where is the anti-Hermitian group generator. Using 



and 



1 211 r <Pp 



z i J (27r)2 p 



(4.12) 



(4.13) 



we can write ( |4.12|) in [Euclidean] momentum space. The term with n potentials is given by 



SA- 



N + 



6n 



tr 



A+{ki) . . . A+{k„ 



where 



-An 



ip + qo)iP + qi) ■■■{? + Qn) 

i 



F{k,,...,k^) (4.14a) 



(4.14b) 



(4.15) 



Using the identities 



cPp d 
(27r)2 dp 



{P + Qi 



(4.16a) 
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d / 1 ' 
dp \p + q ^ 



ttS'-'^^ {p + q) 



(4.16b) 



we can evaluate F as 



F{ki,k2,...,kn) = J2 



(Qo - qt){.qi - Qi) ■ ■ ■ {.Qt-i - - qi) 



In - qi 



(4.17) 



We can now continue to Minkowski space by A; — » A; ■ Q, k k ■ Q'. Since the A's have e~*'^^ 
factors, the retardation condition requires k^ ~^ + 'i^i- We thus define F in Minkowski 
space as the expression ( [4.17|) with 



qi = ^{kj ■Q + 

i 



3) ' 



The expression for the current, which follows, is therefore written as Wl 

d'^ki d'^kn 



(27r)4 (27r)^ 



^(J:>^y^jl^ (k) 



(4.18) 



(4.19a) 



N + 



+ 



(2^: 



Qtt 

.n+l 



dQ 



tr 



tr 



A_{k^) + A+{k,] 



A+{k{)...A+{kn) 



F{ki, ...,kr, 



6n,l (4.19b) 



The equations of motion for the potential with the hard thermal loop corrections is 

{D^F'^na = J'a (4.20) 

with from (4.19). (Evidently we have not displayed gauge-fixing and ghost terms.) Equa- 
tion ( [4.20| ) describes Landau damping of fields in the quark-gluon plasma as well as Debye 
screening and propagation of plasma waves in a gauge covariant way. [Of course, as men- 
tioned in the Introduction, these equations still do not include the soft momentum loop- 
integrations which have to be done in terms of effective vertices and propagators defined by 
equations (4.19), {WM} -\ 

The two-point function [i.e. t?, = 1) as defined by (4.19) agrees with our previous consid- 
erations. The three-point function {i.e. n = 2) as given again by (4.19) is 



/ ]\f \ F ^ f 



tr 



A+{k2)A+{k^) 



X 



{k2 -Qks-Q' -k2-Q'ks-Q) 
{k2-Q + ie){k3 -Q + ie) [{k2 + ks) ■ Q + ie] 



(4.21) 
{Q Q') ] 
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For quark loops we do not have to consider T*-products or contact terms like since 
j^l,a _ iqJ^^^T°■qJ_ Equation ( [4 .81 ) holds as it is written and we can calculate the 3-point 
function (i.e. two A's) by evaluating the retarded commutators. One can verify that this 
does indeed give equation (4.21) with the ie's as indicated. 

Our emphasis has been on the operator equations of motion and on how retarded functions 
arise in such a description. There are of course many other real time formalisms for thermal 
field theory. For example, in the Keldysh approach [§] one considers time to run along a 
contour, which passes from — oo to oo, and then folds back and returns from oo to — oo. We 
have checked that the evolution equations for field configurations in this formalism lead to 
the retarded functions as in our discussion. 

Prescriptions for obtaining the imaginary parts of Green's functions in terms of 'cutting 



rules' have been discussed in Ref. 1112 
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